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Abstract 

The study of short-duration gamma-ray bursts (GRBs) has undergone a revolution in 
recent years thanks to the discovery of the first afterglows and host galaxies in May 2005. 
In this review we summarize our current knowledge of the galactic and sub-galactic environ- 
ments of short GRBs, and the implications for the progenitor population. The most crucial 
results are: (i) some short GRBs occur in elliptical galaxies; (ii) the majority of short GRBs 
occur in star forming galaxies; (iii) the star forming hosts of short GRBs are distinct from the 
host galaxies of long GRBs in terms of star formation rates, luminosities, and metallicities, 
and instead appear to be drawn from the general field galaxy population; (iv) the physical 
offsets of short GRBs relative to their host galaxy centers are significantly larger than for 
long GRBs; (v) the observed offset distribution agrees well with predictions for the locations 
of NS-NS binary mergers; and (vi) unlike long GRBs, which tend to occur in the brightest 
regions of their hosts, the environments of short GRBs generally under-represent the light 
distribution of their host galaxies. Taken together, these observations suggest that short 
GRB progenitors have a wide age distribution and generally track stellar mass rather than 
star formation activity. These results are fully consistent with NS-NS binary mergers, but 
partial contribution from prompt or delayed magnetar formation is also consistent with the 
data. 



1 Introduction 

One of the most fundamental questions in the study of astrophysical transient phenomena is 
the identity of their underlying progenitor systems. The connection between transients and 
progenitors provides insight into the energy source and the explosion/eruption mechanism, events 
rates and their evolution over cosmic time, and a census of the birth and destruction rates of 
compact objects (white dwarfs, neutron stars, and black holes). Ideally, a unique association 
between transient events and their progenitors may come from the identification of individual 
progenitor systems in pre-explosion observations. For example, t he progenitor of sup ernova 



SN 1987A in the Large Magellanic Cloud was a blue supergiant (e.g.. lGilmozzi et al.lll987l ). while 



the progenito rs of several type IIP supernovae (SNe) have been identified as red supergiants (e.g., 



Smarttll2009i ). Unfortunately, for a wide range of cosmic explosions, the low rate of nearby events 
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and /or the intrinsic faintness of the likely progenitor systems preclude a direct identification at 
the present. 

In such cases we may still gain deep insight into the nature of the progenitors through 
statistical studies of the galactic and local environments of the explosions. For example, past 
studies of supernova environments have demonstrated that Type la and Type Ib/Ic/II events 
arise from distinct progenitor systems since the former are located in a wide range of galaxy types, 
while the latter occur only in sta r forming galaxies, pointing to a direct link with massive stars 



(e.g., van den Bergh et al. 20051 ). In a similar vein, long-duration gamma-ray bursts (GRBs) 



were initially linked with massive stars through their e xclusive association with star forming 
galaxies (e.g., Bloom et al. 19981 : Djorgovski et al. 1998 : Fruchter et al. 19991 ). This connection 
was subsequently confirmed by the p resence of accompanying type Ic core-collapse SNe (e.g., 
Hjorth et all 1200.4 IStanek et al.1 120031 ). 

In this review we focus on the study of the galactic and sub-galactic environments of short- 
duration GRBs (T90 ^ 2 s), and discuss the implications for their still unknown progenitors. 
This study is in its nascent stages. Short GRB afterglows were only discovered in 2005 and the 
discovery rate of new events is modest (~ 10 per year). Still, the sample is now large enough 
that we can begin to examine the environments of short GRBs, particularly in comparison with 
long GRB hosts and field galaxy samples, as well as with th eoretical expectations for popular 



proge nitor models such as NS-NS and NS-BH binary mergers (jEichler et al.lll989l : iNarayan et al 
1992). In £}2]we discuss the identification of short GRB host galaxies and measurements of their 



redshifts. The galactic-scale properties (luminosities, star formation rates, and metallicities) are 
described in Finally, in ^]we focus on the locations of short GRBs within their host galaxies, 
partly in comparison to long GRBs and NS-NS merger predictions. The basic results of these 
studies are that the progenitors of short GRBs are distinct from the massive star progenitors of 
long GRBs, that they are related to an older stellar population, and that short GRB progenitors 
reside in average environments within a representative sample of galaxies, likely selected by mass 
rather than by star formation. 

The material presented in this review d raws heavily f r om thr e e recent pape rs on the proper 



Bereer et all (|2007bl l. [Be7geTI (j2009h . and lFong et al 



ties of short GRBs and their host galaxies 

(|2009h - and we refer the reader to these papers for additional details. 



2 Host Galaxy Identifications and Redshifts 

The discovery of afterglow emission from short GRBs starting in May 2005 led to the first 
identifications of their host galaxies and hence to redshift measurements. The first short burst 
to be localized to a positional accuracy of a few arcseco nds, GRB 050509b, appeared to coincide 
with the outskirts of an elliptical galaxy at z = 0.226 (IGehrels et al.ll2005l : iHjorth et al.ll2005al : 
Bloom et al1l2006l ). However, the statistical significance of the association was only ~ 10 3 , and 
the afterglow error circle contained several fainter galaxies possibly at higher redshift. Only two 
months later, the afterglow of GRB 050709 was localized to a sub-arc second position coincident 
with the outskirts of an irregular star forming galaxy at z = 0.161 (|Fox et al.1 [20051 ^ . Despite 
the on-going star formation activity within the host galaxy, the burst was not accompanied 
by a supernqya exp losion, indicating that the progenitor was not likely to be a massive star 
(jHiorth et al.ll2005bh . However, due to the presence of active star formation an association with 
a young progenitor system such as a magnetar could not be excluded. 

The discovery of X-ray, optical, and radio afterglow emission from GRB 050724 finally estab- 
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Figure 1: Discovery of the near- infrared afterglow of GRB 070724. In each row we display the afterglow 
imag e, template image, a nd residual image. The afterglow coincides with the disk of an edge-on galaxy. 
From lBerger et all (|2009f ). 



lished a direct link with an old stellar population (jBerger et al.ll2005T ). The burst was localized 
to an elliptical galaxy at z = 0.257 with no evidence for on-going star formation (<^ 0.05 Mq 
yr -1 ) and a stellar population age of ^ 1 Gyr. The absence of star formation activity and an 
associated supernova demonstrated a connection with an old stellar population. 

The combination of low redshifts (z ~ 0.1 — 0.3) and apparent dominance of early-type host 
galaxies in the early sa mple of short GR Bs led to ini tial claims of a particularly o ld progenitor 
popul ation: r > 4 Gyr iWakar et alfe pOfiT) , r £ 7 Gyr rtzheng fc Ramirez-R,uiz!l2007r ) . r ~ several 
Gyr ( Gal- Yam et al. 20081 ) . Indeed, some a uthors have noted a possible inconsistency with the 
delay time distribution of NS-NS binaries kakar et al.l lioOrV although pop ulation synthesis 



mode ls of NS-NS binary formation and mergers have led to opposite claims (jBelczynski et al 
20060 . Clearly, the sample of short GRBs with afterglow detections that was available when 



these various claims were published was very small (4 events: GRBs 050509b, 050709, 050724, 
and 051221). 

Fortunately, the continued detection of short GRBs, primarily by Swift, and a community- 
wide concerted effort to detect and study their afterglows led to a substantial increase in the 
sample (e.g., Figures CD and [2]) , and a re-evaluation of the host galaxy demographics and redshift 
distribution. For com pre hensive s umma ries of the host galaxy properties we refer the reader to 
Berger et al] (|2007bh and|BergerJ (|2009h . The first comprehensive study o f this sample , and th e 
implications for the redshift distribution of short GRBs, was presented in lBerger et al.l (|2007bl ) . 
Using optical follow-up observations of nine short GRBs with afterglows available as of the end 
of 2006 we found that eight are likely associated with faint galaxies, R ~ 23 — 26.5 mag. By 
comparison to the early host galaxies (with R ~ 17 — 22 mag and z <^ 0.5), as well as the hosts of 
long GRBs and large field galaxy samples, we demonstrated that these new host galaxies likely 
reside at z ~ 1; see Figure [3l Indeed, spectroscopic redshifts for the four brightest hosts led to 
z ~ 0.4 — 1.1. The current redshift distribution of short GRBs (in comparison to long GRBs) is 
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Figure 2: Ground-based images from Magellan and Gemini of several short GRB hosts. All images are 
20" on a side, with the exception of GRB 060502b which is twice as large. The large circles mark the 
XRT error regions, while smaller circles mar k the position s of the optical afterglows (when available). 
Arrows mark the positions of the hosts. From lBerger et aT] (|2007bf ). 



shown in Figured! 

These observations established for the first time that 1/3 — 2/3 of all short GRBs originate 
at z ^ 0.7, and that some bursts produce 10 50 — 10 52 erg in their prompt emission, at least 
two orders of magnitude larger than the low redshift short bursts. Most importantly, with this 
new high redshift sample, we found tighter constraints on the progenitor age distribution than 
previously possible. Viable models include a wide log-nor mal distribution with r* ~ 4 — 8 Gyr, 
or a power law distribution, P(r) oc r n , with — 1 <^ n <J (jBerger et al.ll2007bl ). 



3 Host Galaxy Luminosities, Metallicities, and Star Formation 
Rates 

The association of some short GRBs with elliptical galaxies demonstrated unambiguously that 
at least some of the progenitors are related to an old stellar population. However, as we discussed 
in the previous section, a substant ial fraction of short G RBs (1/3 — 2/3) reside at higher redshifts 
than previously suspected, z ~ 1 (jBerger et al.ll2007bh . and spectroscopic observations indicate 
that most of these galaxies are undergoing active star formation. Indeed, in the sample of 
short GRBs localized to better than a few arcseconds about 50% reside in star forming galaxies 
compared to only ~ 10% in elliptical galaxies; the remaining ~ 40% are currently unclassified 
due to their faintness, a lack of obvious spectroscopic features, or the absence of sensitive follow- 
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Figure 3: Host galaxy R magnitudes as a function of redshift for short GRBs (solid black circles, arrows), 
long GRBs (gray squar es), and galaxies in the HST/ ACS Early Release Observation fields VV 29 (UGC 
10214) and NGC 4676 (Be nitez et al.l[2003 ). The upward pointing t riangles indicate the median redshift 
of a galaxy sample complete to the appropriate magnitude limit (jCoe et al.ll2006l ). For R J> 26 mag 
appropriate for our sample the median redshift is about 1.1. The bottom panels show the redshift 
distributions of galaxies in two magnitude bins from spectros copic (23 < R < 25 mag; ICowie et al.l 
12004 IWirth et all I2004T ) and photometric (25 < R < 27 mag; ICoe et al.l [20061) redshift surveys. The 
clear magnit ude-redshift relation for short GRB hosts suggests that the faint host galaxies are located at 
z ~ 1. From lBerger et al.l (|2007bD . 



up observations. This result raises the question of whether some short GRBs are related to star 
formation activity rather than an old stellar population, and if so, whether the star formation 
properties are similar to those in long GRB host galaxies. The answer will shed light on the 
diversity of short GRB progenitors. 

As of December 2007, the sample of host galaxies with spectroscopic observations was com- 
prised of 6 systems with sub-arcsecond afterglow positions, and an equal number based on only 
Swift/XKT positions (2 — 5 arcsec). The distribution of absolut e rest-frame B-band magnitudes 
(Mb) for these 12 host galaxies ranges from about 0.1 to 1.5 (|Bergerll2009h . with the exc eption 



of the possible elliptical host galaxy of GRB 050509b with Lb ~ 5 ( Bloom et al. 2006). 



The star form ation rates inferred from the [OII]A3727 line using the standard conversion 
( Kennicutt 19981 ). indicate values of 0.2 — 6 M yr _1 , while the elliptical hosts have upper limits 
of ^ 0.1 M yr _1 . Combined with the absolute magnitudes, we find specific star formation rates 
of SFR/Lb « 1 - 10 M yr" 1 L~ l for the star forming hosts (and £ 0.03 M yr" 1 L~ l for the 
elliptical hosts). The SSFR values as a function of redshift are shown in Figure 03 

Finally, for five of the twelve host galaxies we have sufficient spectral information to 
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Figure 4: The redshift distribution of long and short GRBs as of late 2009. 



measure the metallicit y. We use the standard metallicity diagnos tics, R23 = (-^[oii]A3727 + 



[OIII]AA4959,5007j 



)/Fnp (|Pagel et alJll97Sl : iKobulnickv Kewlevll2004h and F [NII ] A6584 /FH a . The 



value of R23 depends on both the metallicity and ionization state of the gas, which we determine 
using the ratio of oxygen lines, 32 = F[oiii]AA4959,5007/F[Oii]A3727- W e note that the R23 dia; 



Kewlev fc Dopital 1200 



nostic is double-valued with low and high metallicity branches (e.g 
This degeneracy can be broken using [Nil] /Ha when these lines are accessible. To facilitate a 
subsequent comparison with field galaxies we use the R23, O32, and [Nil] /Ha calibrations of 
Kobulnickv &: Kewlevl (|2004l ). The typical uncertainty inherent in the calibrations is about 0.1 
dex. 

Adopting the solar metallicity from Asplund et al. ( 20051 ) , 12 + log(0/H) = 8.66 we find 
12 + log(0/H) ~ 8.6 for the upper R23 branch and ~ 8.0 — 8.5 for the lower branch for the host 
of GRB 061006. For the host of GRB 070724 we find 12 + log(0/H) rj 8.9 for the upper branch, 
and ~ 7.6—8.1 for the lower branch. We find a similar range of values for the host of GRB 061210, 
but the ratio -F]nh] /-^Ha ~ 0.2, indicates 12 + log(0/H) ^ 8.6, thereby breaking the degeneracy 
and leading to the upper branch solution, 12 + log(0/H ) 8.9. For the host of GRB 051221a 
we use the line fluxes provided in ISoderherg etaD (SB), and derive similar values to those for 
the host of GRB 070724. Finally, for the host galaxy of GRB 050709 we lack a measurement 
of the [Oil] emission line, and we thus rely on [Nil] /Ha to infer 12 + log(0/H) ~ 8.5. The 
dominant source of uncertainty in this measurement is the unknown value of O32, but using a 
spread of a full order of magnitude results in a metallicity uncertainty of 0.2 dex. For the hosts 
with double-valued metallicities (GRBs 051221a, 061006, and 070724) we follow the conclusion 
for field galaxies of similar luminosities and redshifts that the appropriate values are those 
for the R93 upper branch (Kobulnickv Sz Kewlev 20041 ) . This conclusions was advocated by 
Kobulnickv fc Kewlevl (|2004l ) based on galaxies in their sample with measurements of both R23 
and [Nil] /Ha. It is similarly supported by our inference for the host galaxy of GRB 061210. 
The metallicities as a function of host luminosity are shown in Figure [6l 
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Figure 5: Specific star formation rates as a function of redshift for the host galaxies of short GRBs 
(black squares), long GRB s (blue circles) and field galaxies from the GOODS-N survey (red crosses; 
iKobulnickv k. Kewlevll2004D . Upper limits for the elliptical hosts of GRBs 050509b and 050724 are also 
shown. The cross-hatched region marks the median and standard deviation for the long GRB host sample. 
The inset shows the cumulative distributions for the three samples. The K-S probability that the short 
and long GRB hosts are drawn from the same distributio n is only 0.3 %, while the strong overlap with 
the field sample leads to a K-S probability of 60%. From iBergerl (|2009f l. 



To place the host galaxies of short GRBs in a broader context we compare their proper- 
ties with those of long GRB hosts and field star forming galaxies from the GOODS-N survey 
(jKobumickv fc Kewlevll2004l ). ^n terms of absolute magnitudes, the long GRB hosts range from 
Mr ^ -15.9 to - 21.9 mag, with a median value of {M B ) rj -19.2 mag ((Lb) ~ 0.2 L 



Berger et al.ll2007ah . Thus, the long GRB hosts extend to lower luminosities than the short 



GRB hosts, with a median value that is about 1.1 mag fainter. A Kolmogorov-Smirnov (K-S) 
test indicates that the probability that the short and long GRB hosts are drawn from the same 
underlying distribution is 0.1. On the other hand, a comparison to the GOODS-N sample reveals 
a similar distribution, and the K-S probability that the short GRB hosts are drawn from the 
field sample is 0.6. 

We reach a similar conclusion based on a comparison of specific star formation rates. For 
long GRB hosts the inferred star formation rates range from about 0.2 to 50 M yr _1 , and 
their specific star formation rates are abou t 3 — 40 M Q yr _1 L~ , with a median value of about 
10 Mq yr _1 L^ 1 ( Christensen et al. 2004] ) . As shown in Figure El the specific star formation 
rates of short GRB hosts are systematically lower than those of long GRB hosts, with a median 
value that is nearly an order of magnitude lower. Indeed, the K-S probability that the short 
and long GRB hosts are drawn from the same underlying distribution is only 3.5 x 10~ 3 . This 
is clearly seen from the cumulative distributions of specific star formation rates for each sample 
(inset of Figure [5]). On the other hand, a comparison to the specific star formation rates of the 
GOODS-N field galaxies reveals excellent agreement (Figure [5]). The K-S probability that the 
short GRB hosts are drawn from the field galaxy distribution is 0.6. Thus, short GRB hosts are 
drawn from the normal population of star forming galaxies at z <^ 1, in contrast to long GRB 
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Figure 6: Metallicity as a function of £?-band absolute magnitude for the host galaxies of short GRBs 
(black squares) and long GRBs (blue circles). The gray bars mark the 14 — 86 percentile range for 
galaxies at z ~ 0.1 from the Sloan Digital Sky Survey llTremonti et al. 2004 ) , while r ed crosses designate 
the same field galaxies at z ~ 0.3 — 1 shown in Figure O (Kobul nickv fc K cwlcv 20Q3). Both field samples 
exhibit a cle ar luminosity-metall icity relation. The long GRB hosts tend to exhibit lower than expected 
metallicities (jStanek et alj 12006). while the hosts of short GRBs have higher metallicit ies by about .6 
dex, are moreover in excellent agreement with the luminosity-metallicity relation. From Berger (2009). 



hosts, which have elev ated specific star formation rates, likely as a r esult of preferentially young 



starburst populations ( Christensen et al. 2004 ; Savaglio et al. 20081 ). 



Finally, the metallicities measured for short GRB hosts are in excelle nt agreement with the 
luminosity-metallicity relat ion for field galaxies at z ~ 0.1 — 1 (Figure [6J Kobulnickv &: Kewlev 
IQjj; iTremonti et~ai1l2004h . The two hosts with M B ps -18 mag have 12 + log(0/H) « 8.6, 
while those with Mb ~ —20 to —21 mag have 12 + log(0/H) « 8.8 — 8.9, following the general 
trend. On the other hand, the short GRB host metallicities are systematically higher than 
those of long GRB hosts, which have been argued to have lower than expected metallicities 



fetaoek et ahl B. The median metallicity of short GRB hosts is about 0.6 dex higher than 



for long GRB hosts, and there is essentially no overlap between the two host populations. 

To conclude, the short GRB host sample is dominated by star forming galaxies, but these 
galaxies have higher luminosities, lower star formation rates and specific star formation rates, 
and higher metallicities than the star forming host galaxies of long GRBs. Instead, the short 
GRB host sample appears to be drawn from the typical field galaxy population represented for 
example by the GOODS survey. These results suggest that while short GRB hosts are mainly 
star forming galaxies, the progenitor population most likely traces stellar mass rather than star 
formation activity. 
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4 Hubble Space Telescope Study of the Galactic and Sub-galactic 
Environments of Short GRBs 



High angular resolution imaging of short GRB host galaxies can provide detailed inf ormation on 



the ho st morphologies and for the first time, the burst sub-galactic environments. In lFong et al 
(2009) we performed the first comprehensive analysis of HST observations of short GRB host 
galaxies. We used a sample of 10 events covering the period of May 2005 to December 2006. 
Of these 10 bursts, seven have been localized to su b-arcsecond prec ision and of those, six are 
robustly associated with host galaxies (for details see &^D B . In two cases the identity 



of the host remains unclear. Illustrative examples of the host images and morphological model 
fits are shown in Figure [71 To study the sub-galactic environments we used the offsets relative 
to the host center and the fractional flux at the GRB position. 
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Figure 7: Top: HST/WFPC2/F814W image of the host galaxy of GRB 051221 with a ha error circle 
representing the afterglow position. Center: Sersic model f it from galfit. Right: Residual image. 
Bottom: Same, but for the host galaxy of GRB 061006. From lFong et all (|2009h . 



4.1 Morphological Analysis 

We modeled the two-dimensional surface brightness distributions of the host galaxies to deter- 
mine their effective radii and morphological properties such as the Sersic n index. We found that 
three hosts (GRBs 050709, 051221a, and 061006) are best modeled with n 1, corresponding 
to an exponential disk profile, while two hosts (GRBs 050509b and 050724) are best modeled 
with n ~ 3 and ~ 5.6, respectively, typical of elliptical galaxies. The final three hosts (GRBs 
051210, 060121, and 060313) are equally well modeled with a wide range of n values, although 
n ~ 1 is preferred. Therefore, of the eight short GRB host galaxies with HST observations only 
two can be robustly classified as elliptical galaxies based on their mo rphology. A similar fraction 
was determined independently from spectroscopic observations (JJJ Berger 20091 ) . 



The morphological analysis also yields values of the galaxy effective radii, r e . We find a 
range of i=a 0.2 — 5.8", corresponding to physical scales^ of about 1.4 — 21 kpc. The smallest 



1 For the faint hosts without a known redshift (GRBs 051210, 060121, 060313, and possibly 061201) we assume 
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effective radius is measured for the host of GRB 060313, while the host of GRB 050509b has the 
largest effective radius. The median value is r e ~ 3.5 kpc. The effective radii as a function of 



n ar e shown in Figure 



from IWainwright et al. 



Also shown are the same r e and n values for the hosts of long GRBs 
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Figure 8: Effective radii for the short GRB hosts with HST observations plotted as a function of their 
Sersi c n values. Also shown are the data for long GRB hosts based on HST observations (jWainwright et all 
120071 ). The hosts of GRBs 050509b and 050724 have n values typical of elliptical galaxies, but the 
remaining hosts have a similar distribution to that of long GRBs (i.e., a median of n ~ 1 ; or an exponential 
disk profile). On the other hand, the hosts of short GRBs are large r by about a facto r of 2 than the hosts 
of short GRBs, in agreement with their higher luminosities. From Hong et al.l (|2009l ). 

Two clear trends emerge from this comparison. F irst, all long GRB host s have n <^ 2.5, 
and the median value for the population is (n) ~ 1.1 (jWainwright et al.ll2007l ). Thus, all long 
GRB hosts are morphologically classified as exponential disks, while 2 of the 8 short GRB hosts 
studied here exhibit de Vaucouleurs elliptical galaxy profiles. However, for the hosts with n <J 2, 
the distributions of n values for both populations appear to be similar. 

Second, short GRB hosts have l arger effective radii, wi th (r e ) ~ 3.5 kpc, compared to 



(r e ) « 1.7 kpc for long GRB hosts (jWainwright et al.1 120071 ). A K-S test indicates that the 



probability that the short and long GRB hosts are drawn from the same underlying distribution 
of host galaxy effective radii is only 0.04. Thus, we conclude with high significance that short 
GRB host galaxies are systematically larger than long GRB hosts. The large r sizes of short GRB 
hosts are expected in the context of the galaxy sizeduminos ity relation (e.g., Freeman 1970l ) and 
the higher luminosity of short GRB hosts ($3j Berger 20091 ) . 

An additional striking difference between the hosts of long and short GRBs is the appar- 
ent dearth of interacting or irregular galaxies in the short GRB sample. Of the eight short 



z = 1 l|Berger et al.ll2007bh . and take advantage of the relative flatness of the angular diameter distance as a 
function of redshift beyond z ~ 0.5. 
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GRB host galaxies only one exhibits an irregular morphology (GRB 050709) and none appear 
to be undergoing mergers. In contrast, the fracti on of long GRB hosts with an irregular or 
merger /inter action morphology is about 30 — 60% ( Wainwright et all 2007 ). The interpretation 
for this high merger /inter action fraction in the long GRB sample is that such galaxies are likely 
undergoing intense star formation activity triggered by the merger /inter action process, and are 
therefore suitable sites for the production of massive stars. The lack of morphological merger 
signatures in the short GRB sample indicates that if any of the hosts have undergone significant 
mergers in the past, the delay time betw een the merger and the production of a short GRB is 
^ 10 9 yr (e.g., Barnes Hernquist 19921 ). 



4.2 Offset Distribution 

Based on the astrometric tie of the HST host observations to ground-based afterglow observa- 
tions, we found that the projected offsets of the bursts relative to their host galaxy centers are 
in the range of ~ 0.12 — 17.7". The corresponding projected physical offsets are about 1 — 64 
kpc, with a median value of about 3 kpc. The largest offsets are measured for GRBs 050509b 
and 051210, but these are based on Siuz/i/XRT positions only, with statistical uncertainties of 
about 12 and 18 kpc, respectively. If we consider only the bursts with sub-arcsecond afterglow 
positions we find that the largest offset is 3.7 kpc (GRB 050709), and that the median offset for 
the 6 bursts is 2.2 kpc. In the case of GRB 061201 the host association remains ambiguous, but 
even for the nearest detected galaxy the offset is about 14.2 kpc. 

To investigate the offset distribution in greater detail we supplement the values measured 
from the ffS T observations with offsets for GRBs 070724, 071227, and 090510 from ground-based 
observations dBerger et al.l l200fll : Eaii et al]l200d ). In the case of GRBs 070 724 and 071227 the 
optical afterglows coin cide with the disks of edge-on spiral galaxies (Figure [TJ Berger et al. 20091 ; 
D'Avanzo et alJl200fll ). The offsets of the three bursts are 4.8, 14.8, and 5.5 kpc, respectively. 
There are 7 additional events with optical afterglow identifications. Of these bursts, tw o 



(070707 and 070714b) coincide with galaxies (jPiranomonte et al.1 liooi ; iGraham et alJl2009l ). 



but their offsets have not been measured by the respective authors. Based on the claimed 
coincidence we conservatively estimate an offset of <^ 0.5", corresponding to ^ 4 kpc. Two 
additional bursts (070809 and 080503) do not have coincident host galaxies to deep limits, but 
the nearest galaxie s are l ocated about 6.5 and 20 kpc from the afterglow positions, respective!}!! 



(jPerlev et alJl200«L boOflh . For the final three bursts (080905, 090305, and 090426) no deep host 



galaxy searches exist in the literature. 

In addition to the bursts with sub-arcsecond positions, several hosts have been identified 
within XRT error ci rcles in follow-up observations (G RBs 060801, 061210, 061217, 070429b, 
070729, and 080123; iBerger et ai1l2007bl : iBergerl l200fll ) . Since the putative hosts are located 
within the error circles , in all of these cas es the offsets are consistent with zero or may be as 
large as ~ 30 kpc (e.g., Berger et al.ll2007bh . For example, the offsets for GRBs 060801, 061210, 
and 070429b are 19 ± 16 kpc, 11 ± 10 kpc, and 40 ± 48 kpc. We use 30 kpc as a typical upper 
limit on the offset for these 6 events. We note that no follow-up observations are available 
in the literature for most short GRBs with X-ray positions from 2008-2009. Finally, about 



2 GRB 070714b is located at z = 0.923, while the redshift of GRB 070707 is not known. Based on the faintness 
of the host, R « 27.3 mag, we assume z = 1 to calculate the physical offset. 

3 GRB 070809 is located 19.6 kpc from a galaxy at z = 0.219, and about 2.3" from a much fainter galaxy, which 
at z > 1 corresponds to 18.4 kpc. No host is detected at the position of GRB 080503 in deep HST observations, 
but a faint galaxy is located about 0.8" away, which at z :> 1 corresponds to 6.5 kpc. 
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Figure 9: Projected physical offsets for short GRBs (black) and long GRBs (gray: [Bloom et al]|2002l) . 
The top panel shows a cumulative distribution, while the bottom panel shows the differential distribution 
taking into account the non-Gaussian errors on the offsets. The arrows in the bottom panel mark the 
median value for each distribution. The median value for short GRBs, rj 5 kpc, is about a factor of 5 
times larger than for long GRBs. The arrows in the top panel exhibit the most robust constraints on the 
offset distribution, taking into account the fraction of short GRBs with only 7-ray positions, as well as 
short GRBs for which hosts have been identified within XRT error circles (thereby providing a typical 
range of ~ — 30 kpc). Also shown in the top panel are predicted offset distributions for NS-NS binary 
mergers in Milky Way type galaxies based on population synthesis models. We find good agreement 
betw een the observed d istribution and models, as well as between the robust constraints and models. 
From lFong etaU (|2009f ). 



1/4 — 1/3 of all short GRBs discovered to date have only been detected in 7-rays with positional 
accuracies of a few arcminutes, thereby precluding a unique host galaxy association and an offset 
measurement. 

The cumulat ive distribution o f projected physical offsets for the GRBs with HST observations 
from our work (jFong et alJ l200flh . supplemented by the bursts with offsets or limits based on 



optical afterglow positions (070707, 070714b, 070724, 070809, 071227, 080503, and 090510) is 
shown in Figured! Also shown is the differential probability distribution, P(5r)d(5r), taking 
into account the non-Gaussian errors on the radial offsets (see discussion in Appendix B of 



Bloom et al.1 12002 ). We find that the median for this sample is about 5 kpc. 



As evident from the discussion above, this is not a complete offset distribution; roughly an 
equal number of short GRBs have only limits or undetermined offsets due to their detection in 
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just the X-rays or 7-rays^l Taking these events into account, our most robust inferences about 
the offset distribution of short GRBs are as follows: 

• At least 25% of all short GRBs have projected physical offsets of 10 kpc. 

• At least 5% of all short GRBs have projected physical offsets of ^ 20 kpc. 

• At least 50% of all short GRBs have projected physical offsets of <J 30 kpc; this value 
includes the upper limits for the hosts identified within XRT error circles. 

These robust constraints are shown in Figured 

We compare the observed distribution and the robust constraints o utlined above with 
predicted distrib u tions for NS-NS binar ies in Milky Way type galaxies ( Bloom et al. 19991 : 



Frver et al.1 Il999l : iBelczvnski et al.1 l200d ). appropriate for the observed luminosities of short 



GRB host galaxies (IBergerl 120091) . We fin d go od agreement between t he observed distribution 



a nd those predicted by Bloom et al. ( 19991 ) and Belczvnski et al. ( 200d ). The offset distribution 
of lFrver et aD (jl999h . with a median of about 7 kpc, predicts larger offsets and therefore provides 
a poorer fit to the observed distribution, which has a median of about 5 kpc. However, all three 
predicted distributions accommodate the offset constraints. In particular, they predict about 
60 - 75% of the offsets to be £ 10 kpc, about 80 - 90% to be £ 30 kpc, and about 10 - 25% 
of the offsets to be ^ 20 kpc. Thus, the projected physical offsets of short GRBs are consistent 
with population synthesis predictions for NS-NS binaries. However, the observations are also 
consistent with partial contribution from other progenitor systems with no expected progenitor 
kicks, such as WD- WD binaries. 

We compare ou r observed short GRB offsets with those of long GRBs from the sample of 
Bloom et al.1 (120021 ) in Figure M The offset distribution of long GRBs has been used to argue 
for a massive star progenitor population, and against NS-NS binaries (jBloom et al.ll2002l ). The 
offset distribution for short GRBs is clearly shifted to larger physical scales. In particular, the 
median offset for the long GRBs is 1.1 kpc, about a factor of 5 times smaller than the median 
value for short GRBs. Similarly, no long GRB offsets are larger than about 7 kpc, whereas at 
least some short GRBs appear to have offsets in excess of 15 kpc. 

In the context of NS-NS binary progenitors, the close similarity in the normalized offset 
distributions can be interpreted to mean that most systems likely remain bound to their hosts 
(rather than ejected into the intergalactic medium), and/or have a relatively short delay time. 
These conclusions are tentative due to the small number of events with host-normalized offsets, 
but they can be further tested with future HST observations. 



4.3 Light Distribution Analysis 

In addition to the offset analysis in the previous section, we study the local environments of 
short GRBs using a comparison of their local brightness to the host light distribution. This 
approach is advantageous because it is independent of galaxy morphology, and does not suffer 
from ambiguity in the definition of the host center (see Fruchter et al. 200d ) . We note that for 
the overall regular morphology of short GRB hosts the definition of the host center is gen erally 
robust, unlike in the case of long GRBs (jFruchter et al.l l200d : IWainwright et al.l 120071 ) . On 
the other hand, this approach has the downside that it requires precise pixel-scale positional 
accuracy. In our sample, this is the case for only 6 short bursts. 



4 We do not consider the bursts that lack host searches since there is no a priori reason that these events (mainly 
from 2008-2009) should have a different offset distribution compared to the existing sample from 2005-2007. 
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Figure 10: Cumulative distribution of fractional flux at the location of short GRBs relative to their host 
light. For each burst we measure the fraction of host light in pixels fainter than the GRB pixel location. 
The shaded area is defined by the results for the two available filters for each short GRB. Also shown 
are data for long G RBs (d ark gray line) an d for core-collapse and Type la SNe (light gray lines) from 
iFruchter et al.l (|2006f) and lKellv et alj (|2008| ). The dashed line marks the expected distribution for objects 
which track their host light distribution. Short GRBs appear to undcr-rep resent their host ligh t, while 
long GRBs tend to be concentrated in the brightest regions of their hosts ([Fruchter et al.ll2006f ). From 
iFong et all (|2009t ). 



The fraction of host light in pixels fainte r than the afterglow pixel brightness for each 
host/filter combination is given in IFong et al. (|2009l ). The cumulative light distribution his- 
togram is shown in Figure [TUJ The shaded histogram represents the range defined by the dual 
filters for 5 of the 6 bursts. We find that the upper bound of the distribution is defined by the blue 
filters, indicating that short GRBs trace the rest-frame optical light of their hosts better than 
the rest-frame ultraviolet. This indicates that short GRB progenitors are likely to be associated 
with a relatively old stellar population, rather than a young and UV bright population. 

The overall distribution has a median value of ~ 0.1 — 0.4 (red); namely, only in about 
one-quarter of the cases, 50% of the host light is in pixels fainter than at the GRB location. 
Thus, the overall distribution of short GRB locations under-represents the host galaxies' light 
distribution. This is also true in comparison to the d istribution for core-collapse SNe, which 
appear to track thei r host light ( Fruchter et al. 20061 ). and even Type la SNe, which have a 
median of about 0.4 ( Kelly et al. 20081 ). Thus, the progenitors of short GRBs appear to be more 
diffusely distributed than Type la SN progenitors. 

An extensi ve analysis of the bri ghtness distribution at the location of long GRBs has been 
carried out by Fruchter et al. ( 20061 ). These authors find that long GRBs are more concentrated 
on the brightest regions of their hosts than expected from the light distribution of each host. 
In particular, they conclude that the probability distribution of GRB positions is roughly pro- 
portional to the surface brightness squared. As can be seen from Figure [T0l short GRBs have a 
significantly more diffuse distribution relative to the host light than long GRBs. In particular, 
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for the latter, the median light fraction is about 0.85 compared to about 0.25 ±0.15 for the short 
GRBs. 



4.4 Implications for the Progenitors 



Our extensive analysis of short GRB host galaxy morphologies and the burst local environments 
has important implications for the progenitor population. We address in particular the popular 
NS-NS merger model, as well as delayed mag netar formation via WD- WD mergers or WD 



accretion- induced collapse (jMetzger et al.ll2008l ). 



Morphology: From the morphological analysis we find continued evidence that the bulk of short 
GRB host gal axies (~ 3/4 ) are late-type galaxies, in agreement with results from spectroscopic 
observations (IBerger 2009). Moreover, as demonstrated by the systematic differences in lumi- 
nosity, star formatio n rates, and metallicities between the star forming hosts of long and short 
GRBs ( Berger l200flh . we find that short GRB hosts are systematically larger than long GRB 
hosts. These results indicate that the progenitors of the two GRB classes select different envi- 
ronments. The higher luminosities, larger sizes, and lower specific star formation rates of short 
GRB hosts suggest that their rate of occurrence is tied to galactic mass rather than to star 
formation activity. This result is in broad agreement with old progenitor populations such as 
NS-NS, NS-BH, or WD- WD binaries, but it indicates that the bulk of short GRB progenitors 
are not young magnetars. This conclusion is also supported by the dearth of galaxy merger 
signatures, which point to delays of ^> 10 9 yr relative to any merger-triggered star formation 
episodes. 

Offsets: The differential offsets measured from the HST observations provide the most precise 
values to date for short GRBs, with a total uncertainty of only ~ 10 — 60 mas, corresponding 
to ~ 30 — 500 pc. We find that none of the offsets are smaller than ~ 1 kpc, while this is the 
median offset for long GRBs. On the other hand, a substantial fraction of the measured offsets 
are only a few kpc. The median offset for the HST observations supplemented by ground-based 
data is about 5 kpc (Figure [9]), about 5 times larger than for long GRBs. 

As discussed above, the observed offset distribution is incomplete. About 1/4 — 1/3 of 
all short GRBs have only 7-ray positions (~ 1 — 3'), and a similar fraction have only XRT 
positions, which generally lead to a range of offsets of ~ — 30 kpc. Taking these limitations 
into account we find that the most robust constraints on the offset distribution are that ^ 25% 
of all short GRBs have offsets of <^ 10 kpc, and that ^ 5% have offsets of ^ 20 kpc. Both the 
observed offset distribution and these constraints are in good agr eement with predictions for the 



offset distribution of NS-NS binaries in Milky Way type galaxies (jBloom et al.lll999l ; iFrver et al 



19991 : iBelczvnski et allliooi ). However, at the present they cannot rule out at least a partial 
contribution from other progenitor systems such as delayed magnetar formation and even young 
magnetar flare. The apparent existence of large offsets in the sample suggests that these latter 
models are not likely to account for all short GRBs. 

In the context of implications fo r the progenitor population, a recent study of short GRB 
physical offsets by Troia et al. ( 20081 ) led these authors to claim that short GRBs with extended 
X-ray emission have systematically smaller offsets, possibly due to a systematic difference in the 
progenitors. Our HST sample includes three short GRBs with strong extended emission (050709, 
050724, and 061006) , and one GRB (060121) with possible extended emission (4.5cr significance; 
Donaghv et"aT1l2006l ). The physical offsets of these bursts are about 3.7, 2.7, 1.3, and 1 kpc, 
respectively, leading to a mean offset of about 2.2 kpc. The physical offsets of the bursts without 
extended emission, but with precise afterglow positions (051221, 060313, and 061201) are 2.0, 
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2.3, and 14.2 or 32.5 kpc, respectively. The two events with no extended emission and with 
XRT positions (050509b and 051210) have offsets of about 54 ± 12 and 28 ±23 kpc, respectively. 
Including the ground-based sample with optical afterglow positions, we find that the bursts 
with apparent extended emission (070714b, 071227, 080513, and 090510; barbier et all 120071 : 



Sakamoto etaD 120071 : Rjkwatta et alll2009l : IPerlev et al1l2009h have offsets of £ 4, 14.8, ~ 20, 



and ~ 5.5 kpc, while the bursts without extended emission (070724 and 070809) have offsets of 
4.8 and ~ 6.5 kpc. Thus, based on the sample of events with sub-arcsecond positions we find 
that 6/8 bursts with extended emission have offsets of ^ 5 kpc and 2/8 have likely offsets of 
~ 15 — 20 kpc. In the sample without extended emission we find that 4/5 have offsets of ^ 6 
kpc and 1/5 has a likely offset of ~ 14 — 32 kpc. Thus, we conclude that there is no significant 
difference in the two offset distributions. 

The inclusion of events with only XRT positions does not change this conclusion. In partic- 
ular, of the subset with no extended emission only GRB 050509b is likely to have a significant 
offset, while GRBs 051210, 060801, and 070429b have offsets (28 ± 23, 19 ± 16, and 40 ± 48 
kpc, respectively) that are consistent with zero. Similarly, GR B 061210 with exten ded emission 



has an offset of 11 ± 10 kpc. An examination of the sample of Troia et al. (20081) reveals that 



their claim that short GRBs without extended emission have systematically larger offsets rests 
on four events in particular: GRBs 050509b, 060502b, 061217, and 061201. As noted above, 
GRBs 050509b and 061201 indeed appear to have substantial offsets, but so do GRBs 071227 
and 080503 with extended emission and offsets of about 15 — 20 kpc. Next, the large offset for 
GRB 060502b relies on its claime d association with an elliptical galaxy 70 ± 16 kpc from the 
XRT position ( Bloom et al.l 12007). Howeve r, the XRT error circle contains additional galaxies 
with negligible offsets (jBerger et al. 2007b ). Finally, we note that the offset for GRB 061217 
i s unr eliab le due to a substan tial discrepancy of about 33" in the XRT positions from I Butler 
(|2007h and lEvans etaD (j2009h . A continued investigation of the difference between short GRBs 
with and without extended emission will greatly benefit from the use of host-normalized offsets. 

Light Distribution: In addition to projected offsets relative to the host center, we find that the 
locations of the short GRBs with HST imaging and sub-arcsecond positions are more diffusely 
distributed relative to their host light than long GRBs. In particular, we find that short GRB 
positions under-represent their host light, even in comparison to core-collapse and Type la SNe. 
This result is likely an upper limit on the brightness of short GRB locations since only the subset 
of events with optical afterglow positions can be studied with this approach. Thus, short GRBs 
arise from a population of events with a more diffuse distribution than massive stars and Type 
la SN progenitors. This result also indicates that the bulk of the progenitors of long and short 
GRBs cannot both be magnetars. 

There are currently 10 known short GRBs with optical afterglows for which HST observations 
will enable a similar analysis. This is twice the number of the current sample, and we can 
therefore make significant progress in understanding the relation of short GRB environments to 
the overall distribution of light in their host galaxies with future observations. 



5 Conclusions 

While the sample of short GRBs with afterglow positions is still significantly smaller than 
the sample of long GRBs, we have began to make significant progress in understanding their 
galactic and sub-galactic environments. The results of spectroscopic and high resolution imaging 
observations point to an association of short GRBs with faint regions of normal star forming and 
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elliptical galaxies. In nearly every respect (star formation rates, metallicities, sizes, offsets, light 
distribution) the environments of short GRBs are distinct from those of long GRBs, indicating 
that they are not related to a young progenitor population. The bulk of the evidence indicates an 
association with an old stellar population, and the observations are fully consistent with NS-NS 
binary mergers. However, a partial contribution from prompt or delayed magnetar formation 
cannot be ruled out at the present. We expect that studies of the ever-growing sample of short 
GRBs similar to the ones presented in this review will eventually constrain the contributions of 
various progenitors models. 
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